Dynamics of global entanglement under decoherence 
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We investigate the dynamics of global entanglement, the Meyer-Wallach measure, under 
decoherence, analytically. We study two important class of multi-partite entangled states, the 
Greenberger-Horne-Zeilinger and the W state. We obtain exact results for various models of 
system-environment interactions (decoherence). Our results shows distinctly different scaling 
behavior for these initially entangled states indicating a relative robustness of the W state, 
consistent with previous studies. 
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I. INTRODUCTION 

It is well-known that the notion of quantum entangle- 
ment is a key concept in quantum mechanics. It is also 
responsible for "strange" non-local behavior of quantum 
systems in marked contrast to classical notions of reality 
and locality [1,2]. Schrodinger refer to it as the "essence 
of quantum mechanics" [1] . Besides its fundamental as- 
pects, entanglement constitutes the central part of new 
modes of information technology, quantum computation 
and quantum communication [3,4,5] and is therefore a 
key ingredient of many information processing protocols. 
Recently, a considerable amount of work has been de- 
voted to characterize, quantify and realize different vari- 
ety of entangled states [6]. By now, bi-partite entangle- 
ment is a relatively well-understood phenomena. How- 
ever, the situation becomes much more complex when 
multi-partite systems are considered [7] . 

On the other hand, entangled states are very fragile 
when they are exposed to environment. The biggest en- 
emy of entanglement is decoherence which is believed to 
be the responsible mechanism for emergence of the clas- 
sical behavior in quantum systems [8] . Since the mainte- 
nance and control of entangled states is essential to real- 
ization of quantum information processing systems, the 
study of deteriorating effect of decoherence in entangled 
states would be of considerable importance from theoret- 
ical as well as experimental point of view [9,10,11]. 

To demonstrate the effects of decoherence on entangled 
states, an appropriate entanglement measure which could 
be capable of monitoring the dynamics of entanglement 
in decoherence processes is needed . However, there are 
no exact measure of entanglement under general condi- 
tions for mixed states. Even for bipartite mixed states, 
apart from the particular case of two-level systems [12], 
the exact solution is missing. There is an approximate 
generalization of concurrence for mixed states which was 
proposed by Mintcrt et al [13] and has bean used to mea- 
sure entanglement in multi-qubit systems. 

Global entanglement (GE), defined by Meyer-Wallach 
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(MW) entanglement measure of pure-state [14] which is 
a monotone[15] , is a very useful measure of entanglement. 
As we will show briefly, GE is a measure of total non-local 
information per particle in a general multi-qubit system. 
Therefore, GE gives an intuitive meaning to multi-qubit 
entanglement as well being an experimentally accessible 
measure [16]. In this paper, we use this measure to mon- 
itor the entanglement dynamics of two seminal multi- 
qubit entangled states, the Greenberger-Horne-Zeilinger 
(GHZ) state and the W state, under different models of 
system-environment interaction. 

To demonstrate the GE dynamics of multi-qubit entan- 
gled state under decoherence, we need to know the gen- 
eralization of GE to mixed state. Unfortunately, there 
is no generalization of the primary definition of the MW 
measure to mixed states, analytically. However, we can 
monitor the GE dynamics for two important classes of 
multi-qubit entangled state, the GHZ and W state, by 
exploiting the relationship between GE and tangles. To 
elucidate this point, we adopt informational approach 
which can give an intuitive meaning to GE. 



II. GLOBAL ENTANGLEMENT, 
INFORMATION, AND TANGLES 

Finite amount of information can be attributed to N- 
qubit pure state which is N bit of information according 
to Brukner-Zeilinger operationally invariant information 
measure [17]. This information can be distributed in lo- 
cal as well as non-local form, which is associated with 
entanglement [18]. This information has a complimentary 
relation: 

-^total — Ilocal ~t~ I'non—local- (1) 

The total information is conserved unless transferred 
to environment through decoherence. The amount of 
information in local form is Iiocai = X^i^i where, 
li = 2Trpf — 1 is the operationally invariant information 
measure of a qubit [17]. Therefore, according to Eq.(l) 

Inon-iocai = Yl,i=i 2(1 " Trp}) which can be distributed 
in different forms of quantum correlations, the tangles. 
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among the system, 

^non— local — 2 ^ ^ Ti-^^i^ ~\~ ~\~ N ^ ] '^ii...ijv? (^) 
ji<i2 ii<j2<---<iiv 

where the first term is referred to as 2-tanglc, the next 
being 3-tangle and the last term the N-tangle of the sys- 
tem. It can easily be seen that the MW measure of GE 
can be written as: 

^ai = ^[2 E ^'1'^ + N Yl ^^---l- 

ii<i2 ii<i2<--.<iN 

Therefore, the MW measure of GE (Egi) is a measure 
of total non-local information per particle (or average of 

tangles per particles ■^). This property of Egi resembles 
the molar-extensive thermodynamic variables (e.g. heat 
capacity Cp). In fact, we expect Egi to have thermo- 
dynamically relevant and experimentally accessible fea- 
tures in multi-qubit (spin) systems [19]. This property of 
GE in contrast to most other multi-partite entanglement 
measures. It is therefore of considerable interest to in- 
vestigate the dynamics of GE under decoherence in the 
large N limit for the generic, experimentally realizable, 
entangled states. We note that Egi does not give detailed 
knowledge of tangles distribution among the system. For 
example, Egi cannot distinguish between entangled states 
which have equal (r) yet different distribution of tangles, 
like \GHZ)n and \EPR)®f. However, Egi can distin- 
guish between GHZ and W state since their distinctly 
different type of tangles leads to different values of (r) . 

We consider the two seminal multi-qubit entangled 
state: \GHZ) = -i=(|00...0) + |ll...l)) and \W) = 

^(|00...01) -I- I00...10) -I- ... -I- I10...00)) as initial states 

which arc known to bear qualitatively different quantum 
correlation. For example, GHZ state has a different tan- 
gle distribution than that of W state. In operational con- 
text they cannot be transformed to each other by local 
operation and classical communication (LOCC) [20]. So 
it is interesting to investigate the entanglement behav- 
ior of this two type of multi-qubit entangled states under 
different models of system-environment interaction. 

In the GHZ state, all the entanglement is contained in 
the N-tangle form of non-local information which can be 
computed by N-concurrence measure (Cjv) [21]. Cjv is 
a generalization of two-qubit measure, concurrence [11], 
for even number of qubits (N). Therefore, for GHZ state 
GE is Egi = tn = C% . This can be easily computed for 
even N as: 

Rn = pjvaf ^p^af ^, (4) 

2"-l 

Cjv = max{Xi — Aj,0}; 

1=2 

where ay is Pauli matrix in y direction and A^s are the 
eigenvalues of the matrix Rn with Ai being the maxi- 
mum. 



In the W state, only two-qubit entanglement, 2-tangles 
(T2), is present with each qubit equally entangled with all 
other qubits. T2 can be computed analytically from the 
above measure. Thus, there are different two- 

qubit entanglement in the W state. GE is therefore given 
by: 

Egi = {N- l)r2, (5) 

for W state. It is important to note that this approach 
allows us to obtain the exact solution of GE dynamics 
under decoherence of different environment models. This 
also shows the great utility of GE as a measure of entan- 
glement. 

III. DECOHERENCE MODEL SYSTEMS 

In order to evolve our chosen states under influence of 
decoherence we use the Lindblad form of master equation 

[22], 

, N 

f = EA.. (6) 

The Lindblad operators, Li, describe the local interac- 
tion of each qubit with environment independent of other 
qubit interaction with the environment. We assume Lj 
is the same form for all qubit, Li = L. For markovian 
process [22] 

^^'^ = E?PJfcp4-{JfcJl,p}], (7) 
k ^ 

where operator Jj, describes the system-environment 
model of interaction with strength 7^. In this paper we 
investigate dissipative, dephasing, and noise processes, 
each with a well-defined Jfc. For the two- level systems 
the operators, J^, are expressed in terms of Pauli matri- 
ces. The solution of Lindblad form of master equation 
for two- level systems are studied in [23]. 

For dissipative environment, Ji = (T_. In this process 
the system interacts with a thermal bath at zero temper- 
ature. This process could be described as spontaneous 
emission of a two-state atom coupled with the vacuum 
modes of the ambient electromagnetic field which leads 
the atom state to the ground state. For dephasing pro- 
cess, Ji = o"+(T_. This is a phase-destroying process 
that does not have a classical counterpart and is there- 
fore intrinsically quantum mechanical. It corresponds to 
a situation where no energy is exchanged with environ- 
ment, that is, the population of energy eigenstates of the 
system do not change with time. Only the phase informa- 
tion which includes quantum correlations is lost. For the 
noisy environment, J\ = and J2 = C7+. Noisy dynam- 
ics are related to another extreme of thermal bath, i.e. 
when temperature is extremely high while the system- 
bath coupling is extremely weak. This process random- 
izes the state of the system which results in a completely 
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mixed state eventually. The noise process has a partic- 
ular interest since its effect is basis independent. That 
is, the noisy operation is invariant under unitary oper- 
ation. All these processes could have different effect on 
the multi-qubit entangled state. But the common feature 
of them is that under the action of each of these environ- 
ments any initial entangled state asymptotically evolves 
to a separable state. 

IV. RESULTS 



The dynamics of GE in dephasing, dissipation and noisy 
environment for the initial W state is illustrated in Figs.l 
and 2. Eq.(10)(Fig.l) shows that the decay rate(a, for 
Egi oc exp(— at) for GE is independent of N for the W 
state in dissipative and dephasing environment as found 
previously using numerical solution for a different mea- 
sure of entanglement [10]. Note, however, that the rate 
of change of GE decreases with increasing N. For noisy 
environment. Fig. 2, we observe a decay to separable state 
after a finite time t^ep which increases linearly with N, 
also consistent with previous studies [10]. 



Our goal is to obtain the time dependence of the den- 
sity matrix, p{t) — e~^*/5(0), of the system in order to 
determine the time evolution of Egi(t) for the initially 
prepared multi-qubit entangled states, i.e. W and GHZ. 
According to the structure of entanglement in W state 
we can deduce the time dependency of the GE from two- 
qubit entanglement, T2. The two-qubit density matrix, 

(0) = jf\'ip+){^J+\ -f 1 00) (00 1, is the same for any 

pair of qubits, ij. Therefore, for the initial W state in 
the dissipative process we have. 
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where p = e is the decoherence parameter. Similarly, 
for dephasing process, one obtains. 
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and, therefore, for the initial W state, GE has the simple 
exact solution. 



iV2 



(10) 



for both dissipative as well as dephasing processes. For 
noisy process, the density matrix is 



{N-2) 
N 

which leads to 
iV- 1 
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maa;{4p^-(l-p^)2(Ar2_(p7V-4p)'')2,0}]^ 

(12) 
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FIG. 1: Egi vs. "ft in an initial W state with dephasing (or 
dissipative) process for — 2, 4, 6, 8 qubits. The number 
of qubits increases from the top curve to bottom. The inset 
shows the decay rate vs. A^. Note that these are simple graphs 
of Eq.(lO) and are only drawn for comparison with other Figs. 
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FIG. 2: Egt vs. 'yt in an initial W state with noisy process for 
A*' = 2, 4, 6, 8 qubits. The number of qubits increases from the 
top curve to bottom. The inset shows the (linear) dependence 
of decay rate on A'^ for up to A^ = 14 qubits. 

For the GHZ state, the density matrix in the dephasing 



4 



process is 

p^{t) = i(|o)«^(or^ + p^io)«^(ir^ + (13) 

which leads simply to Egi = e~^^*'. The GHZ density 
matrix for the dissipative process is 

p^{t) = \{p^ |o)«^(l|«^ + pf |l)«^(0|«^ (14) 
1 

+ 51 A(z)|gig2...gw)(<7ig2...gAr|), 

(ji,...,giv=0 

where 

N 

■Z' = ^ ; 2i = — — - , 
i=i 

and 

(Z) = {qiq2...qN\Z\qiq2...qN), 

and 

A(z)=b<^>(l-p)^-<^>+0<^>]. 
The GHZ density matrix for the noisy process is 

pj,{t) = i(p^|0)®^(l|«^ + p^|l)«^(0|«^ (15) 
1 

+ X! A(2)|gig2---gAr)('?ig2---gAr|), 

(ji,...,(jjv=0 

Consequently, our results (Egi), for the GHZ state in 
dephasing, dissipative and noisy environment are shown 
respectively in Figs. 3, 4, and 5 for various N. As shown 
in the corresponding insets, the decay rates increases lin- 
early with system size N for all three processes. Also, 
the rate of change of Egi increases with system size as 
well. All these results are consistent with previous stud- 
ies using different methods than ours[9,10,ll]. For the 
dynamics of GE in GHZ state, although we have the ex- 
act solution only for even number of qubit, the behavior 
of GE under decoherence for the odd number can be in- 
ferred from the simplicity and symmetry of our results. 
For example, our results for Egi in dephasing process 
holds for any number of qubit in the initial GHZ state. 

V. CONCLUSION 

In conclusion, in this work we have given an intuitive 
informational meaning of MW measure of global entan- 
glement. Based on the relationship between global en- 
tanglement and tangles, which constitute the non-local 
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FIG. 3: Egi vs. 'yt in an initial GHZ state with dephasing 
process for A'^ = 2, 4, 6, 8 qubits. The number of qubits in- 
creases from the top curve to bottom. The inset shows the 
(linear) dependence of decay rate vs. A'^ up to A'^ = 10 qubits. 
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FIG. 4: Egi vs. 'yt in an initial GHZ state with dissipation 
process for A'' = 2, 4, 6, 8 qubits. The number of qubits in- 
creases from the top curve to bottom. The inset shows the 
(linear) dependence of decay rate on iV up to Af = 10 qubits. 



form of the information, we identify the exact solution of 
its dynamics under different system-environment models 
for the two qualitatively different multi-qubit entangled 
states, the GHZ and W states. In all the cases consid- 
ered, we obtain an exponential decay of entanglement as 
a function of time. For the W state, the results show 
that the lifetime of the GE is independent of the number 
of the qubits in dephasing and dissipative processes and 
the lifetime linearly decreases with N in a noisy process. 
While for the GHZ state, the lifetime of GE decreases 
linearly with N. Our results indicate that the quantum 
correlations in W state are more robust to decoherence 
effects than that of the GHZ state. 

The authors kindly acknowledge the support of Shiraz 
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FIG. 5: Egi vs. "/t in an initial GHZ state with noisy process 
for A'^ = 2, 4, 6, 8 qubits. The number of qubits increases 
from the top curve to bottom. The inset shows the (linear) 
dependence of decay rate vs. A'^ up to A*' = 10 qubits. 
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